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Abstract: Aeolian processes have been studied extensively at low elevations, but have been relatively little 
studied at high elevations. Aeolian sediments are widely distributed in the Yarlung Zangbo River basin, 
southern Tibetan Plateau, which is characterized by low pressure and low temperature. Here, we 
comprehensively analyzed the wind regime using data since 1980 from 11 meteorological stations in the 
study area, and examined the interaction between the near-surface wind and aeolian environment. The wind 
environment exhibited significant spatial and temporal variation, and mean wind speed has generally 
decreased on both annual and seasonal bases since 1980, at an average of 0.181 m/(s-10a). This decrease 
resulted from the reduced contribution of maximum wind speed, and depended strongly on variations of 
the frequency of sand-driving winds. The drift potential and related parameters also showed obvious spatial 
and temporal variation, with similar driving forces for the wind environment. The strength of the wind 
regime affected the formation and development of the aeolian geomorphological pattern, but with variation 
caused by local topography and sediment sources. The drift potential and resultant drift direction were two 
key parameters, as they quantify the dynamic conditions and depositional orientation of the aeolian 
sediments. Wind affected the spatial variation in sediment grain size, but the source material and complex 
topographic effects on the near-surface wind were the underlying causes for the grain size distribution of 
aeolian sands. These results will support efforts to control aeolian desertification in the basin and improve 
our understanding of aeolian processes in high-elevation environments. 
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1 Introduction 


The Tibetan Plateau is one of the Earth's most important and unique geographical units, and is 
experiencing serious desertification due to its cold and semi-arid to arid climate, with less than 500 
mm of annual precipitation (Han et al., 2009; Dong et al., 2017). Due to its high elevation, at an 
average of more than 4500 m above sea level (a.s.l.), its low temperatures, and the existence of a 
cryosphere, aeolian sediments have developed extensively in the plateau's dry basins, floodplains, 
river terraces, lakeshores, and mountain slopes (Dong, 2017). As a result, the plateau is an 
important dust source region that has profound impacts on the ecology, atmosphere, and marine 
environments of East Asia and even the whole northern hemisphere (Fang et al., 2004; Xu et al., 
2018). 

The Yarlung Zangbo River passes through a typical area of aeolian sedimentary deposits in the 
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southern Tibetan Plateau. Due to favorable environmental conditions, including the abundant sand 
source, strong winds, and presence of many locations suitable for sediment accumulation (Li et al., 
1999), aeolian sediments, including loess deposits and dune sands, have been deposited in and 
around the basin, especially in the basin's wide valleys (Here, we defined wide valleys as valleys 
in which the valley floor, including the river terraces, was at least 3-4 km wide, versus a width of 
less than 3 km in most of the river's valley). Under the influence of both natural factors and human 
activities, the basin's aeolian sandy land has gradually expanded, at an annual growth rate of 764.7 
hm? from 1975 to 2008 (Shen et al., 2012). During the winter and spring dry seasons, the riverbed 
and floodplain sediments in the wide valleys are entrained and transported by the wind, resulting 
in strong aeolian activity, and cause widespread damage to the plateau's ecological environment, 
while endangering air transport and the production and daily life of local residents. For example, 
Zhongba County, in the upper reaches of the river, has experienced three relocations since the 
county's establishment in 1960 due to the impact of aeolian activity (Sun et al., 2010). Several 
studies have shown that although human disturbance since the 1950s has significantly accelerated 
the development of sandy desertification in several parts of the Tibetan Plateau, the extent of these 
disturbances has been small (due to the plateau's small and concentrated population) and their 
effects have therefore been overwhelmed by the impacts of climate change, which is the main 
driving force behind evolution of the ecological environment (Dong et al., 1999; Zou et al., 2002; 
Wang et al., 2016; Dong, 2017; Dong et al., 2017). 

Wind is the most important force responsible for the development of desertification, and plays 
an important role in the erosion, transport, and accumulation of surface sediments (Dong et al., 
2011). It is also a critical climate factor that affects aeolian activity and desertification. Many 
scholars have used remote sensing to monitor the evolution of sandy desertification in the Yarlung 
Zangbo River basin (Li et al., 2004; Sun et al., 2010; Yuan et al., 2010; Li et al., 2011; Shen et al., 
2012). They have discussed the evolution of the sandy land and its possible climate driving factors, 
such as temperature and precipitation, but our literature review found no in-depth analysis of the 
basin's wind regimes. Only variations of the mean wind speed have been discussed. Although the 
mean value can reflect the general trend for a wind environment, it cannot provide enough detail 
to understand aeolian movement because the relationship between wind speed and sand transport 
intensity is nonlinear (Dong et al., 2011). Sand drift potential is a good index to assess the intensity 
of aeolian activity, as this variable represents the wind's potential ability to transport sediment; as 
a result, it has been widely used in aeolian research. Cui et al. (2017) discussed the spatiotemporal 
variation of the wind energy environment in the dune fields of northern China, and provided a 
reference for regional desertification control. Yang et al. (2019) showed that drift potential was 
closely related to the migration rate of sand dunes, and was an important parameter that reflected 
the intensity of aeolian activity. It is therefore essential to provide an integrated analysis of the sand 
drift potential in the basin. However, we found insufficient details in the literature to allow such an 
analysis. 

To provide key missing information, we combined field investigations with wind data from 11 
weather stations in and around the Yarlung Zangbo River basin from 1980 to 2015 (the most recent 
data available), and analyzed the spatial and temporal variation of the wind environment and sand 
drift potential. Based on the results of this analysis, we discussed the relationships between the 
wind regimes and the basin's aeolian geomorphological pattern and sand transport patterns. This 
will expand the data available to aeolian researchers, and assist the prevention and control of 
aeolian hazards in the basin. More importantly, these results provide insights into the formation and 
development of surface dust (i.e., loess and sand dune) in the southern Tibetan Plateau, and 
strengthen our understanding of aeolian processes in high-elevation arid and semi-arid areas. 


2 Physiographic setting 


The Yarlung Zangbo River basin (28°59'—30°34' N, 82°32'—94°28' E) stretches across most of the 
southern Tibetan Plateau from west to east, and the average elevation of the riverbed is close to 
4000 m a.s.l. (Fig. 1). The valley is an area of the plateau with deep negative terrain, with the river 
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cutting to a depth of 0 to 150 m below the surrounding land in the upper reaches and up to 2000 m 
in the lower reaches (Li et al., 1997; Shen et al., 2012). Wide and narrow valleys alternate, and 
loess deposits and sand dunes have developed throughout the basin (Yang et al., 2020). The loess 
has mostly accumulated since the last glacial maximum, with development in different periods 
since the late glacial (Ling et al., 2020). The sand dunes are relatively young, and mostly formed 
hundreds or even thousands of years ago (Li et al., 1997). Typical aeolian landforms such as 
barchan dunes and barchanoid chains have developed, mainly on the valley floors (Li et al., 1999). 
Climbing dunes, sand sheets, and loess deposits have mostly developed on the valley slopes (Sun 
et al., 2007). These deposits have mainly developed in the four wide valleys shown in Figure 1: 
from west wo east, these are the Maquanhe, Shigatse, Shannan, and Mainling wide valleys. The 
deposits cover a total area of 2737 km’, of which 55.5% represent mobile and semi-fixed states 
(Shen et al., 2012). The aeolian activity is strong and the surface dust can be blown to a height of 
more than 1000 m above the valley floor, especially in the spring and winter (Fig. 2). 
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Fig. 1 Tibetan Plateau and its major atmospheric systems (a), and the distribution of the main sandy land, the 
location of sample sites (numbered S1 to S6 from west to east) and meteorological stations in and near the Yarlung 
Zangbo River basin (b). The data from site S6, in the Mainling wide valley, were obtained from Zhou et al. (2011). 
Station names: PR, Purang; TR, Tingri; LZ, Lhaze; SGS, Shigatse; GT, Gyantse; NM, Nyemo; NGZ, Nagarze; LS, 
Lhasa; ZT, Zetang; NC, Nyingchi; BM, Bome. 


The basin lies within a transitional region between an arid to semi-arid climate and a humid 
climate, with precipitation ranging from 251 mm in the west to 580 mm in the east, and a mean 
annual potential evaporation ranging from 2293 mm in the west to 2734 mm in the east (Liu et al., 
2019). The mean annual temperature ranges from 3.9°C to 9.2°C, with mean monthly temperatures 
in west ranging from —7.3°C in January to 14.8°C in July, and mean monthly temperatures in the 
east ranging from 1.1°C in January to 17.4°C in July (based on daily temperature data from 
meteorological stations in and near the basin from 2001 to 2015). The soil types include alpine 
meadow soils, alpine steppe soils, mountain shrubby steppe soils, and grassland soils. The 
vegetation is dominated by Sophora moorcroftiana, Stipa purpurea, and Orinus thorodii. Due to 


202101.00057v1 


chinaXiv 


ChinaX ivé (ERAT! 


JOURNAL OF ARID LAND 


the dry and cold climate, the fragile ecological environment, and the low vegetation cover, wind 
erosion is intensified and aeolian activity is frequent (Liu et al., 2019). 


(b) (c) 


Fig. 2 Visible evidence of aeolian activity in the Yarlung Zangbo River basin. The blue arrows indicate the 
location of the river, and the white arrows indicate the direction of sand transport. E, east. Photos were taken by Dr. 
YANG Junhuai from 15 to 28 November 2019. 


3 Data and methods 


3.1 Data 


In this paper, we selected 11 meteorological stations in and near the Yarlung Zangbo River basin 
(Fig. 1), and obtained daily-scale data from the climate database provided by the Chinese National 
Meteorological Information Center (CMIC, http://data.cma.cn/en). To ensure consistency of the 
statistical period among the stations, we used data from 1 January 1980 to 31 December 2015. The 
raw data included the daily mean wind speed, daily maximum wind speed, and the corresponding 
wind directions. The daily mean wind speed was the average of the values measured at 02:00, 
08:00, 14:00, and 20:00 local time (the data were recorded in Beijing standard time). The maximum 
wind speed based on the 10-min average speed. All data were collected at a height of 10 m above 
the ground surface, as recommended by the World Meteorological Organization. 

We used the mean and maximum wind speeds to characterize variations in the wind environment. 
Given the importance of the threshold velocity for entrainment of surface particles (see Section 3.3 
for details), we selected three parameters for each wind speed: the total wind speed, the sand- 
driving wind speed, and the frequency of sand-driving wind (wind frequency). The total wind speed 
represents the time-weighted sum of all measured wind speeds, and was used to analyze the overall 
wind environment. The sand-driving wind speed represents a wind speed greater than the threshold 
speed, and that can therefore entrain sediment. The wind frequency is the proportion of the total 
study period (number of days) in which a sand-driving wind was recorded. To better analyze the 
spatial variation of the wind environment, we used data from the Purang station to describe the 
wind regime in the Maquanhe wide valley; data from the Tingri, Lhaze, Shigatse, and Gyantse 
stations to describe the wind regimes in the Shigatse wide valley; data from the Nyemo, Nagarze, 
Lhasa, and Zetang stations to describe the wind regimes in the Shanan wide valley; and data from 
the Nyingchi and Bome stations to describe the wind regimes in the Mainling wide valley. The 
Purang, Tingri, Lhaze, Gyantse, Nyemo, Nagarze, and Bomi were missing data on the maximum 
wind speed and direction from 1980 to 2006; this data was available for the other four stations. To 
facilitate comparison and analysis despite the missing data, we focused on the sand drift potential 
from 2007 to 2015. 


3.2 Quality control and homogeneity test 
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These 11 stations belong to the national standard meteorological stations, and these data have been 
passed by the strict quality control. To ensure quality of the data, CMIC carried out a series of 
quality inspections for the data from 1980 to 2010, including climatic limit and allowable value 
check, extremum check, internal consistency check, temporal consistency check, spatial 
consistency check and manual inspection and correction, and adopted strict three-level (station 
level-provincial level-national level) quality control for data from 2011 to 2015. After that, the 
CMIC confirmed that the invalid data was less than 1% and the accuracy of the data was close to 
100%, which suggests the quality of data in a certain extent. 

However, long-term climatic data may be subject to inhomogeneities resulting from artificial 
changes in instruments and observing environment, which will result in uncertainty in the 
assessment of climate trends (Azorin-Molina et al., 2014; Cui et al., 2017). To further test the 
credibility of wind speed data series, we adopted the standard normal homogeneity test (SNHT) 
method (Alexandersson, 1986), which has been widely used in detecting inhomogeneities and 
breakpoints in climate series (Liu, 2000; Wu et al., 2018a). First, we established reference series 
according to the criterion of He et al. (2012). Second, a reference function Q; was established, and 
then computing a new standardized series Z;. Finally, we constructed testing statistic Tk for 
homogeneity test. The maximum value of Tk was defined as To. The threshold of To is 8.10 (n=36, 
P<0.05) according to critical value given in Liu (2000). If To>8.10, the breakpoint is included in 
the data set. Otherwise, the data set is homogeneous. More detailed descriptions of the SNHT 
method can be found in previous studies (Liu, 2000; Wu et al., 2018a). Our results showed that the 
To values range from 4.97 to 7.44, and are less than the threshold of 8.10. Therefore, we considered 
that these wind speed data used in this study are homogeneous and credible. 


3.3 Determination of the threshold wind velocity 


The threshold wind velocity is the critical speed required for sand transport to begin, and is essential 
for predicting sand transport, dust release, and desertification (Li et al., 2019). Some researchers used 
a value of 5 m/s, mostly for sand seas at a much lower elevation than in the present study, as the 
threshold wind velocity in the basin (Liu et al., 2019). However, we assumed that this velocity is 
inconsistent with the actual threshold wind velocity. Bagnold (1941) suggested that the threshold 
wind velocity is dominated by the effects of the grain size of sand particles. However, field 
experiments using a portable wind tunnel showed that the air density is another important factor that 
affects sand particle movement (Han et al., 2014; Han et al., 2015). All else being equal, the wind's 
ability to transport sand will decrease at low air density, and the value of Bagnold's empirical 
parameter A (a proportionality coefficient) will vary with air density (Eq. 3; Han et al., 2015). 

To determine the threshold wind velocity, we collected 39 surface samples (to a depth of 2 cm) 
from typical barchan dunes in the four wide valleys (selected 1 to 2 dunes in each valley), and 
measured the grain size distribution using Malvern Mastersizer 2000 (Malvern Instruments Ltd., UK) 
at the Key Laboratory of Western China's Environmental Systems, Ministry of Education, Lanzhou 
University. The sediments for a barchan dune were collected from the windward toe, the lower and 
upper parts of the windward slope, the crest, the upper and lower parts of the leeward slope, and the 
leeward toe of barchan dunes and interdune areas. We used the average value of grain sizes of the 
samples in each wide valley for analysis. In addition, we obtained estimates of the mean air density 
in the wide valleys (Shi et al., 2019). We calculated the threshold wind velocity (U+) as follows: 


U, =5.75U" xlog Z / Z,. (1) 
U` = Ax (((p, - Pa) / p,)*9 xd), (2) 
A=0.05 + 0.039, , (3) 


where U” is the threshold friction velocity (m/s); Z is the height (10 m) at which we measured the 
wind speed; Zo is the aerodynamic roughness length, which is approximately equal to a grain size 
of 1/30 m; ps is the density of sand grains (2650 kg/m°); pa is the density of air (kg/m°); g is the 
acceleration due to gravity (9.80 m/s”); d is the mean grain size (mm) at each sample site; and A is 
a proportionality coefficient that varies with air density. Because the basin's geographical 
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environment is complex, we calculated the actual threshold wind velocity separately for each of 
the four wide valleys (Table 1). 


Table 1 Mean elevation, air density, coefficient A, and mean grain size of dune sands in the wide valleys of the 
Yarlung Zangbo River basin, and corresponding threshold velocity for sand entrainment 


Wide valley name 


Parameter 
Maquanhe Shigatse Shannan Mainling 
Mean elevation (m) 5000 4000 3700 3000 
Air density (kg/m*) 0.74 0.82 0.85 0.91 
Coefficient A 0.072 0.075 0.076 0.077 
Mean grain size (mm) 0.29 (n=14) 0.24 (n=14) 0.25 (n=7) 0.13 (n=41)" 
Threshold shear velocity (m/s) 8.0 7.2 7.3 5.4 


Note: n represents the number of sand samples; *, data from Zhou et al. (2011). To better reflect the effect of wind on sand particles, we 
excluded samples from the interdune area when calculating the mean grain size. A value is proportionality constants in Equation 3 that 
depend on the air density, and were obtained from Han et al. (2015). 


3.4 Calculation of drift potential 


To better reflect the details of aeolian sand movement, we calculated the following parameters of 
the wind energy environment using the methods of Fryberger and Dean (1979) (Table 2): the drift 
potential (DP), the resultant drift potential (RDP), the directional variability (RDP/DP), and the 
resultant drift direction (RDD). 


DP =U*(U-U,)t, (4) 


where U is the measured maximum wind speed (knots) at 10 m above the ground surface; U; is the 
threshold wind velocity (knots) at 10 m (the wind speed required to entrain sediment; we used the 
value in Table 1); and t is the proportion of the recorded time during which the wind velocity was 
greater than the threshold velocity. More detailed descriptions of the method can be found in Zhang 
et al. (2015). In addition, we defined the seasons as follows: spring, from March to May; summer, 
from June to August; autumn, from September to November; and winter, from December to 
February of the following year. 

To provide more insights into the spatial distribution of DP, we performed kriging interpolation 
to generate contour lines that covered the entire basin. To do so, we used a digital map of the basin 
provide by Google Earth, then interpolated between the meteorological stations using the Spatial 
Analyst Tools feature of version 10.2 of the ArcMap GIS software (www.esri.com). DP represents 
the total annual wind power, and was used to describe the potential capacity of sand transport from 
each wind direction. The value of the RDP/DP close to 1 indicates a single dominant drift direction, 
whereas the value close to 0 indicates multiple significant drift directions. RDD represents the 
resultant direction of the sand transport. These parameters are traditionally presented in a wind rose 
diagram. Although this approach makes a number of interpretative simplifications and assumptions, 
previous studies of the wind energy environments of global sand seas have demonstrated its value 
(Zhang et al., 2018a, b). 


Table 2 Classification of wind energy environments proposed by Fryberger and Dean (1979) 


DP (VU) Wind energy environment RDP/DP Directional variability Wind direction category 
<200 Low <0.3 High Complex or obtuse bimodal 

200-400 Intermediate 0.3-0.8 Intermediate Obtuse or acute bimodal 
>400 High >0.8 Low Wide or narrow unimodal 


Note: DP, drift potential; VU, vector units; RDP, the resultant drift potential. 


4 Results and analysis 


4.1 Overall characteristics of the wind environment 


Table 3 summarizes the mean annual values of the three wind parameters for mean and maximum 
wind speeds during the study period in the Yarlung Zangbo River basin. The results show 
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substantial spatial and temporal differences in the basin's wind environment. 


Table 3 Mean value of wind parameters for mean wind speed (Umean, 1980-2015) and maximum wind speed 
(Umax, 2007—2015) in the wide valleys of the Yarlung Zangbo River basin 


Wide valley name 


Parameter Time Maquanhe Shigatse Shannan Mainling 
V, V, F V, V, F V, V, F V, V, F 
Annual 3.33 90° 0.5 21> 83> 0.4 2.0° 81> 0.2 164 59° 0.1 
Spring 3.3 2.9 0.2 3.0 7.9 0.8 2.4 7.8 0.3 1.9 6.1 0.3 
Unean Summer 3.7 2.9 0.2 1.7 - - 1.7 2.4 0.0 1.4 - - 
Autumn 3.2 8.4 0.3 1.7 2.7 0.02 16 - - 1.4 1.9 0.03 
Winter 3.0 9.5 1.1 32 8.5 0.9 2.1 8.3 0.6 1.7 5.8 0.2 
Annual 7.07 9.0 30.0 7.0 9.2? 42.8 5.6 84° 168 48° 64% 35.0 
Spring 7.4 9.1 37.9 8.1 9.4 63.0 6.3 8.5 26.8 5.3 6.5 49.0 
Umax Summer 7.9 8.9 46.9 6.4 85 311 5.4 8.3 114 4.5 6.3 27.7 
Autumn 6.6 8.9 16.2 5.9 8.6 26.4 4.9 8.2 7.7 4.4 6.2 24.1 
Winter 6.1 9.2 18.9 7.4 9.7 50.0 5.6 8.4 20.9 4.8 6.4 39.0 


Note: V,, total wind speed (the time-weighted sum of all measured wind speeds, m/s); Vs. sand-driving wind speed (m/s); F, the frequency 
of wind with a speed greater than the entrainment threshold (% of days in the study period); -, the wind was not strong enough to entrain 
sediment. Differences between wide valleys for the annual values were compared by means of one-way ANOVA using the SPSS software 
(https://www.ibm.com/analytics). When the ANOVA result was significant, we compared pairs of wide valleys using least-significant- 
difference tests. Values of V, and V; labeled with different letters (a, b, c and d) differed significantly between valleys (ANOVA followed 
by least-significant-difference tests, P<0.05). 


Based on the annual data, the total wind speed, sand-driving wind speed, and its frequency for 
the mean and maximum wind speeds differed significantly among the valleys, and generally 
decreased from west to east. The values for the mean wind speed in the Maquanhe wide valley were 
significantly higher than that in the Shigatse wide valley. However, there was no significant 
difference in the total wind speed values for the maximum wind speed between the Maquanhe and 
Shigatse wide valleys. This is because the mean wind speed can represent the overall mean 
variations in the wind environment, whereas the maximum wind speed represents only a part of 
those overall variations. 

Except for the Maquanhe wide valley, the values of the three parameters were high in spring and 
winter, and were generally the highest in the spring, and were generally low in summer and autumn. 
In the Maquanhe wide valley, the total wind speed for the mean and maximum wind speeds was 
the highest in summer and the lowest in winter, whereas the sand-driving wind speed was the lowest 
in summer and the highest in winter. This suggests that the effect of the wind on sediments will be 
concealed by the measurements of total wind speed, causing researchers to miss some details of 
aeolian sand movement. In addition, the three wind parameters showed relatively consistent 
seasonal variations, which suggests that the sand-driving wind speed and its frequency may be the 
two important factors affecting the strength of the wind environment, but the extent of their 
influence is not yet clear. 


4.2 Wind speed trends 


The overall mean wind speed showed a statistically significant decrease from 1980 to 2015 
(P<0.01), but the magnitude of the decrease differed among the wide valleys (Fig. 3). The 
decreasing trend for the mean wind speed was the greatest in the Maquanhe wide valley (0.25 
m/(s-10a)), and was weakest in the Mainling wide valley (0.13 m/(s-10a)). Except for the Maquanhe 
wide valley, the temporal variation in mean wind speed can be divided into two periods. During the 
first period from 1980 to 2004, the mean wind speed was decreasing. During the second period 
from 2005 to 2015, the mean wind speed began increasing, at a rate of 0.7 m/(s-10a) in the Shigatse 
wide valley (P<0.01), 0.5 m/(s-10a) in the Shannan wide valley (P<0.05), and 0.3 m/(s-10a) in the 
Mainling wide valley. The seasonal trends for the mean wind speed also showed a significant 
decreasing and seasonal differences (Fig. 4a—d). Generally, the maximum rate of decrease occurred 
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in the spring at an average of 0.25 m/(s-10a), and the minimum occurred in the summer at an 
average of 0.13 m/(s-10a). However, the variation trends and seasonal differences for the maximum 
wind speed were not obvious (Fig. 4e—h), possibly because the maximum wind speed only reflects 
part of the overall information about wind speed changes. 
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Fig. 3 Temporal variations in three parameters for the mean wind speed (1980-2015) (a-d) and the maximum 
wind speeds (2007-2015) (e—h) in the wide valleys of the Yarlung Zangbo River basin 
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Fig. 4 Seasonal variations in the mean wind speed (1980-2015) and the maximum wind speed (2007-2015) in 
the wide valleys of the Yarlung Zangbo River basin. k is the slope of the regression and represents the long-term 
trends for the wind speed; R? represents the regression goodness of fit; *™* is significant at P<0.01; and * is significant 
at P<0.05. 


The overall trends for the mean wind speed were roughly consistent with the trends for the 
maximum wind speed from 2005 to 2015 (Figs. 3 and 4), which suggests that the maximum wind 
speed significantly affected the wind environment. To discuss the contribution of the maximum 
wind speed to the overall variations in the wind environment, we used the empirical formula of Cui 
et al. (2017): 


Ciar = Tsx / AT nean ? (5) 
where Cmax is the contribution of the maximum wind speed (%); and Tmax and Tmean represent the 
trends (i.e., regression slopes) for maximum and mean wind speeds, respectively (m/(s-a)). 4Tmean 


is based on the daily mean wind speed calculated from daily four measurement times. It should be 
noted that the equation is valid based on a similar mapping relation between the mean and 
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maximum wind speeds obtained from the wind data records, and may not be valid for datasets with 
other characteristics (Cui et al., 2017). This contribution differed among the stations; it averaged 
62.6% and ranged from 47.6% to 94.5% (Table 4). That is, the variation of the maximum wind 
speed accounted for 62.6% of the total variation of the wind environment on average. 


Table 4 Contribution of the maximum wind speed to the variation in the wind environment (Cmax) 


Meteorological station Cmax (%) 
Annual Spring Summer Autumn Winter 
Shigatse 47.6 94.5 51.5 65.2 48.9 
Lhasa 90.6 65.8 50.8 49.0 50.8 
Zetang 59.5 60.3 53.7 60.6 66.9 
Nyingchi 66.1 48.8 58.5 94.5 69.1 


Note: Data analyzed for the above four meteorological station were for the period 1980-2015. 


Figure 3 also shows that the variation of the total wind speed usually followed the variation of 
the sand-driving wind speed and its frequency, but it was not clear which of these factors played 
the major role. To quantify their effects on the wind environment, we performed multiple linear 
regression to compare the magnitude of the regression coefficients for the sand-driving wind speed 
and wind frequency for the mean and maximum wind speeds. The regression coefficients for the 
sand-driving wind speed ranged from —0.104 to 0.139 (P>0.05), whereas the coefficients for the 
wind frequency ranged from 0.557 to 0.989 (P<0.05). Therefore, the trend for the wind 
environment was only significantly affected by variations of the frequency of the sand-driving wind 
speed. 


4.3 DP 


Although the variations of the mean and maximum wind speeds reflect the overall characteristics 
of the Yarlung Zangbo River basin's wind environment, they provide no detailed information about 
sand transport and wind direction. Based on the contribution of the maximum wind speed, we 
calculated DP and related parameters. 

4.3.1 Spatial variation in DP 

The wind energy environment differed among the four wide valleys. In the Maquanhe wide valley, 
DP was between 200 and 400, so this region had an intermediate wind energy environment. In the 
Shigatse wide valley, DP was far greater than 400, so this region had a high wind energy 
environment. In the Shannan and Mainling wide valleys, DP was less than 200, so these regions 
had a low wind energy environment. Annual RDP/DP (directional variability) was the lowest in the 
Shigatse wide valley, with a wide unimodal wind regime, whereas it was intermediate in the other 
wide valleys, with obtuse or acute bimodal wind regimes, which also explains the differences in 
DP. RDD indicated that the direction of sand transport ranged from north to southeast, and moved 
clockwise (closer to the southeast) from the west end of the basin to the east (Fig. 5). 
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Fig. 5 Wind energy environment in the four typical wide valleys of the Yarlung Zangbo River basin. DP, drift 
potential; RDP, resultant drift potential; RDP/DP, directional variability; RDD, resultant drift direction. 
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To provide further insights into the wind regimes, we analyzed the wind regimes of the 11 
meteorological stations (Fig. 6). The Tingri and Lhaze stations had a high wind energy environment, 
especially at Tingri station (DP=2163 VU). Purang, Gyantse, and Nagarze stations had an 
intermediate wind energy environment. DP was less than 200 at the other stations, indicating that 
they had a low wind energy environment. Although the wind energy environment differed among 
the stations, DP generally decreased from west to east. This is more visible based on kriging 
interpolation of the DP values (Fig. 6). The annual directional variability (RDP/DP) was mostly 
intermediate, but values ranged from 0.1 to 0.9. Nyemo station had the highest directional 
variability, whereas Tingri and Bome stations had low directional variability. RDD values indicated 
that 72.7% of the sand transport direction was towards the northeast. The direction of sand transport 
at Bome station was southeast, versus north-northwest at Nyemo and west at Nyingchi. RDP/DP 
and RDD showed that the direction of near-surface wind in the basin and its variability were 
relatively spatially consistent, but retained clear local characteristics. 
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Fig. 6 Wind energy environment at 11 meteorological stations in and near the Yarlung Zangbo River basin, and 
contour intervals for the DP calculated by kriging interpolation 


4.3.2 Temporal variation in DP 


The data from the long time series at four meteorological stations showed that except for RDP at 
Nyingchi, DP and RDP in the Yarlung Zangbo River basin decreased significantly from 1980 to 
2015 following an exponential equation (Fig. 7). However, there were no statistically significant 
trends at any of the four stations for RDP/DP and RDD. The variation in RDP/DP was relatively 
small at Shigatse and Zetang stations, with annual values of 0.71 (+0.10) and 0.57 (+0.13), 
respectively. RDP/DP at Lhasa increased significantly (R’=0.20, P<0.01), indicating the drift 
direction changed from multiple directions to a single direction. RDP/DP at Nyingchi first 
increased (R*=0.50, P<0.01) and then decreased (R?=0.36, P<0.05), indicating that the sand drift 
direction changed obviously. The variation in RDD was smallest at Shigatse, with a value of 242° 
(+11°). The variation in RDD was moderate at Lhasa and Zetang, with values of 263° (+30°) and 
207° (+31°), respectively. RDD at Nyingchi decreased significantly (R?=0.53, P<0.01), and the 
standard deviation was largest, with a value of +93°, indicating complex sand drift directions in the 
lower reaches of the river during the study period. 
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Fig. 7 Temporal variations in the wind environment at the four meteorological stations with a long statistical 
period (1980-2015). The dashed lines in this figure represent the boundaries for the wind energy regimes defined 
in Table 2. 


We found that the variation in these parameters was similar (Fig. 7). For example, the variations 
in DP and RDP at Lhasa station were more strongly related to those at Zetang station (R°=0.60 and 
0.23, P<0.01, respectively) than to those at Nyingchi station (R?=0.12, P<0.05; R*=0.05, P>0.05, 
respectively). The variation in RDP/DP at Lhasa station was weakly but significantly related to that 
at Nyingchi station (R?=0.17, P<0.05). This suggests that there may be spatial consistency and 
differences in the changes of wind energy environment in the basin, which also supported by the 
variation in the wind regimes in the four wide valleys from 2007 to 2015. For example, there was 
a significant positive correlation between the variations in DP in the Shigatse and Shannan wide 
valleys (R?=0.48, P<0.05), and a weak or even negative correlation among the variation in DP and 
RDP between the Mainling wide valley and the other three wide valleys (Fig. 8a and b). A strong 
positive correlation in RDP/DP existed between the Shannan and Mainling wide valleys (R?=0.59, 
P<0.05; Fig. 8c). The relationships in RDD among the wide valleys were weak and not significant 
(Fig. 8d), which may be caused by the complex topography (discussed in Section 5.2). Overall, the 
characteristics of the parameters showed that the wind energy environment responded to similar 
driving forces in the upper and middle reaches of the Yarlung Zangbo River, especially in the 
middle reaches, but was obviously different in the lower reaches. 

We also observed obvious seasonal variation in the wind energy environment in the basin (Fig. 9). 
DP and RDP were higher in spring and winter than in summer and autumn, except for the Maquanhe 
wide valley (which had high DP in summer rather than winter). RDP/DP was generally lowest in 
winter (i.e., a single dominant drift direction), and highest in summer (i.e., multiple significant drift 
directions), which agrees with the seasonal variation of DP. The seasonal variation in the RDD also 
differed among the wide valleys. In the Maquanhe wide valley, RDD was northeast in spring, and 
rotated counterclockwise towards north-northwest by winter. In the Shigatse and Mainling wide 
valleys, RDD was relatively constant in all seasons, but was east-northeast and east-southeast, 
respectively. In the Shannan wide valley, RDD was south-southeast in summer, but was east- 
northeast in the other seasons. 


202101.00057v1 


chinaXiv 


ChinaXivA ERAT! 


YANG Junhuai et al.: Near-surface wind environment in the Yarlung Zangbo River basin... 


R*=0.87, P<0.01 


120 fo 100 ©) 
1 80 
y= -3.91x+7895.5 5 60 y=-1.84x+3719.1 
20. 


R*=0.83, P<0.01 


0. 0 
2006 2008 2010 2012 2014 2016 2006 2008 2010 2012 2014 2016 


0 
2008 2010 2012 2014 2016 2006 2008 
a Maquanhe e Shigatse a Shannan + Mainling 


2010 2012 2014 2016 


Fig.8 Temporal variations in the wind environment at the four wide valleys (2007—2015). The dashed lines in (c) 
and (d) represent the boundaries for the wind energy regimes defined in Table 2. 
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Fig.9 Seasonal variations in the wind environments of the four wide valleys (2007-2015) 


5 Discussion 


5.1 Variation in the wind environment 


The Yarlung Zangbo River basin is one of the strongest source regions in China for long-distance 
dust transmission, and this dust may have a serious impact both in China and on the global climate 
(Fang et al., 2004; Han et al., 2008). The surface dust can be raised high into the atmosphere even 


JOURNAL OF ARID LAND 


by weak winds and then transported over thousands of kilometers by global wind systems or jet 
streams (Swet et al., 2019). Therefore, surface dust and its dynamic evolution in the basin have 
attracted the attention of many aeolian researchers. However, detailed analysis of the wind 
environment was lacking until the present study. 

The evolution of the mean wind speed shows that the wind environment weakened both annually 
and in each season from 1980 to 2015, at an average rate of 0.181 m/(s-10a). We also found that 
the mean and maximum wind speeds in the middle and lower reaches of the Yarlung Zangbo River 
have increased at different rates since 2005 (Figs. 3 and 4). In fact, this phenomenon that the 
reduction of wind speed from 1980 to 2004 and its recovery after 2005 has been reported in the 
previous studies at global and regional scales. Previous studies suggested that the surface wind 
speed in most parts of the world has declined from mid-1950s until 2000 (Xu et al., 2006; Vautard 
et al., 2010; Wu et al., 2018b; Zhang and Wang, 2020). The most significant decreases have 
occurred in Central Asia and North America, at an average rate of 0.11 m/(s-10a), followed by 
Europe, East Asia and South Asia (0.08 m/(s-10a)), and finally Australia (Wu et al., 2018b). The 
changes in the surface wind speed may be mainly caused by the changes in atmospheric circulation 
(Wu et al., 2018b), land use and cover change (Wu et al., 2016; Zha et al., 2016; Zha and Wu, 2017), 
surface roughness (Vautard et al., 2010), and human activity (Xu et al., 2006). The declining rate 
in this study is higher than the value of 0.140 m/(s-10a) that McVicar et al. (2012) reported for the 
global terrestrial average wind speed, but consistent with the results of Guo et al. (2011) in China 
(0.180 m/(s-10a)), Cui et al. (2017) in the dune fields of northern China (0.193 m/(s-10a)), and Jiang 
et al. (2015) in the northeastern Tibetan Plateau (0.190 m/(s-10a)). This shows that although there 
are regional differences, the wind speeds in China's arid and semi-arid areas have decreased at a 
similar rate. In addition, a reversal of wind speed after 2000 has also been observed at global scale 
(Kim and Paik, 2015; Zeng et al., 2019; Zha et al., 2019; Zhang and Wang, 2020), but the reasons 
for this reversal were complex. This may be partially explained by changes in human activities. 
The population is concentrated in the middle and lower reaches of the river, and increased at a rate 
of 41,000 people per year from 2005 to 2015 (based on data from Tibet's 2018 statistical yearbook). 
More than 80% of local inhabitants live by grazing, farming, and gathering fuel by cutting trees 
(Shen et al., 2012). Decreasing vegetation cover decreases surface roughness (Vautard et al., 2010), 
thereby increasing the near-surface wind speed. We also found the contribution of the maximum 
wind speed to the mean wind speed ranged between 47.6% and 94.5%, which indicates that the 
maximum wind speed played an important role in determining the strength of the wind environment. 
This was consistent with the results of Vautard et al. (2010) and Guo et al. (2011), who found that 
strong winds played a more important role than light winds in defining the wind environment. 

The spatial and temporal characteristics of wind regimes in the different regions in the basin 
revealed that there was both inter-regional consistency and differences in the wind environment. 
The spatial variation in wind regimes indicated that RDP/DP and RDD were relatively consistent 
in the basin, but retained clear local characteristics. The temporal variation in wind regimes 
indicated that the wind environment responded to similar driving forces in the river's upper and 
middle reaches, but was different in the lower reaches. Our results are strongly supported by a 
mathematical simulation of the near-surface flow field in the basin (Liet al., 1999), which indicated 
that differences in the topography changed how the near-surface flow field developed in complex 
ways. Specifically, a westerly wind prevailed in the upper and middle reaches of the river (west of 
the Shannan wide valley), but in the river's middle reaches (near the Shannan wide valley), the 
wind was deflected counterclockwise, which paralleled the changing orientation of the river valley 
in this part of the basin. However, east of the Jiacha mountains (Fig. 1), the valley winds met a 
descending east wind produced by the east Tibet secondary high pressure and downstream back 
flow, and the Indian monsoon can also reach the valley by traveling along the "great gorge" shown 
in Figure 1a (Yang et al., 1989). The interactions among these winds create a complex local flow 
field. 

The temporal variation of the mean wind speed and DP showed that the overall strength of the 
wind environment decreased significantly, but that this decrease did not parallel the evolution of 
the area of aeolian sandy land in the basin. Shen et al. (2012) and Liu et al. (2019) showed that the 
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areas of aeolian sandy lands in the basin (1975 to 2008) and in the four wide valleys (1990 to 2014) 
increased, but that the magnitudes of the increase decreased gradually over time. This may be 
because the weakening wind was less able to entrain surface sediments, thereby allowing the 
sediments to accumulate instead of being transported downwind. Although weakening of the wind 
environment may limit the expansion of aeolian sandy land, the wind remained strong enough to 
provide a dynamic force that has shaped the basin's development since 1975. Previous research 
showed that the development of aeolian sandy land since 1975 was more closely related to human 
activities, increasing air temperature, and decreasing precipitation than to the wind speed (Shen et 
al., 2012). In general, the basin's arid and windy climate created favorable conditions for aeolian 
erosion, whereas human activities (e.g., grazing on natural rangeland) have greatly accelerated the 
development of aeolian sandy land since 1975. 


5.2 Wind regimes and aeolian geomorphological patterns 


DP and its related parameters were related to the valley's aeolian geomorphological patterns. We 
found a moderately strong and statistically significant relationship between DP and the area of 
aeolian sandy land (Fig. 10), which indicates that DP was an important factor that controlled the 
development and spatial distribution of aeolian sand. However, we also found that DP was largest 
in the Shigatse wide valley, but the area of aeolian sandy land was not correspondingly large. This 
suggests that the spatial distribution of aeolian sandy land was at least partially caused by other 
factors. When we excluded the Shigatse data from the regression, we obtained a stronger 
relationship (R°=0.92, P<0.01; Fig. 10). This difference may be because a long and narrow part of 
the valley lies upwind of Shigatse, and the valley changes direction sharply towards the north. This 
topography would have affected wind speed, and thus DP. The sediment supply is also important 
for the development of aeolian sediments (Zhang et al., 2015; Dong, 2017; Zhang et al., 2018b). 
Dong (2017) suggested that available sand sources controlled the characteristics of the sandy land 
on the Tibetan Plateau, such as the distribution and dynamic changes of aeolian sandy land and the 
aeolian landforms. The braided streams that are common in the river valley provide a rich sediment 
supply and an ideal place for the development of aeolian sediments. The river near Shigatse is a 
heavily braided stream 8 to 10 km wide, whereas the braided river near Maquanhe can reach 25 km 
wide, which also explains the results in Figure 10. 
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Fig. 10 The relationship between the DP and the area of aeolian sandy land. The data for the areas in the four 
wide valleys was obtained from Liu et al. (2019). 


RDD strongly affected the sediment distribution patterns. Based on the RDD data for the basin, 
72.7% of the direction of sand transport was towards the northeast (Fig. 6), which was consistent 
with the spatial distribution of dunes in the basin. Aeolian sands developed most often on the slopes 
of the valley, especially along the northern bank of the river (Dong, 2017). Under the action of a 
southwest wind, climbing dunes and sand sheets can remain on the valley's northern slopes at close 
to the angle of repose, and in some cases climbed as far as 6 km along the slope, placing them 
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among the highest-elevation aeolian sands in the world (Dong et al., 2017). The spatial distribution 
of aeolian landforms was also strongly supported by the seasonal variation in RDD (Fig. 9). To the 
west of the Jiacha Mountains, aeolian landforms were mainly distributed on the river's northern 
bank, whereas to the east, they were mainly distributed on the southern bank. 

To further analyze the effect of RDD on sand dunes in the river valley, we analyzed the 
relationship between RDD and dune orientation. We obtained the orientation data by analyzing 
Google Earth images. The local RDD was strongly and significantly (R?=0.99, P<0.01) related to 
the dune orientation, but with the dunes oriented at a mean angle of 2° to 17° to RDD (Fig. 11a). 
However, at some locations, RDD was perpendicular and even opposite to the dune orientation, 
which is not surprising given the basin’s vast expanse and complex topography. The movement of 
the near-surface airflow is easily changed by variations in the local topography, thereby resulting 
in variation in dune orientation (Fig. 11b). Although the overall characteristics of aeolian sand 
transport are shaped by the calculated wind regimes in the basin, our results differ from those 
obtained for lower elevation sand seas in northern China because the complex topography 
prevented us from obtaining a single overall description of sand movement. 
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Fig. 11 Relationship between the RDD and dune orientation (dashed line represents y=x) (a). The dune orientation 
and migration direction (orange arrows) in the Mainling wide valley shown in Figure 1 (center coordinates: 
29°08'43.06”, 93°45'56.28"; Image credit: Google Earth) (b). 


5.3 Wind regimes and sand transport 


The grain size characteristics of the sediments not only reflect the effects of wind on the transport 
and sorting of source materials, but also reflect how land surface characteristics such as topography 
and vegetation affect aeolian sand transport (Zhang and Dong, 2015). The proportion of the 
particles in each grain size category showed that the contents of very fine and fine sands generally 
increased slightly from west to east, but that the contents of medium and coarse sands generally 
decreased (Table 5). That is, the mean grain size in the basin decreased gradually from the upper 
reaches to the lower reaches. The grain size distribution was more strongly related to the sand- 
driving wind speed (R?=0.28 to 0.73) for mean wind speed than to DP (R?=0.007 to 0.32). This 
may be because the mean wind speed (from 1980 to 2015) can better represent the overall 
characteristics of the wind environment, whereas DP is calculated based on the maximum wind 
speed (from 2007 to 2015) and only reflects part of the information on the wind environment. With 
increasing sand-driving wind speed, the proportions of very fine and fine sands decreased, whereas 
the proportions of medium and coarse sands increased (Fig. 12), which agrees with the results of 
Yang et al. (2019) in a series of wind tunnel experiments. Dong and Li (1998) suggested that the 
silt and clay particles (<0.05 mm) are difficult for the wind to erode because of their cohesion and 
because they are more likely to be protected against the wind by larger particles in a deposit that 
has a high proportion of large particles, and this may explain that there was no significant 
correlation between the silt and clay content and the sand-driving wind speed. Our results indicated 
that a strong sand-driving wind will blow away the very fine and fine surface sediments and leave 
medium and coarse sands at the surface; therefore, the grain size of barchan dunes in the valley 
exhibited a certain spatial variation. 
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Table 5 Grain size characteristics of barchan dunes in typical dune fields of the Yarlung Zangbo River basin 


Percentage (%) in each size category (mm) 


Sample Mean size Sorting : 
site  Siltandclay Very fine Fine Medium Coarse (mm) (0) Skewness Kurtosis 
(<0.063) (0.063-0.125) (0.125-0.25) (0.25-0.50) (0.50-2.00) 
S1 6.39 14.41 27.64 35.67 15.88 0.30 1.01 0.21 1.29 
S2 6.74 4.81 42.05 42.22 4.18 0.25 0.90 0.25 1.70 
S3 16.66 13.57 40.02 27.38 2.37 0.19 1.74 0.54 2.19 
S4 4.74 12.61 35.15 35.60 11.90 0.28 0.95 0.10 1.04 
S5 8.36 10.44 38.59 36.43 6.18 0.25 1.20 0.32 1.87 
S6 10.16 30.00 56.47 3.16 0.23 0.13 0.92 0.35 2.07 


Note: The calculation method for the grain size parameters and the related grading standard is provided by Folk and Ward (1957). The 
six sample sites (S1—S6; Fig. 1) run from west to east. The data for the Mainling wide valley were obtained from Zhou et al. (2011). 
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Fig. 12 Relationships between the percentage of different grain size fractions and the sand-driving wind speed 


The sediments in the basin were mostly poorly sorted (59%) and moderately well sorted (33%), 
which indicates that the time and distance of sand transport were relatively short. Of the sediments, 
49% were positively skewed, and 41% were very positively skewed, which indicates that the 
sediments were mainly composed of coarse components. In addition, 56% were very leptokurtic, 
21% were leptokurtic, and 21% were mesokurtic, indicating that the sediment distribution was 
concentrated, and that at least some of the particles directly entered the aeolian environment without 
any change. Several factors are responsible for these characteristics. First, the sand dunes mainly 
developed near or between mountains (such as in the Lhaze area), resulting in short transport 
distances and limits the sorting of aeolian sediments by the wind. Second, previous studies showed 
that the large quantities of coarse fluvial sands that are exposed during the river's dry season 
provided the main material source for sand dunes near the river (Li et al., 1997; Li et al., 1999), 
which agrees well with our description of the sediment deposition environments based on Sahu's 
(1964) hypothesis (Table 6). The barchan dunes in the study area mainly formed in a fluvial 
environment (53.9%) and a shallow marine environment (25.7%). Based on the similarity between 
shallow marine sediments and lacustrine sediments, the shallow marine environment could be 
redefined as a lacustrine environment. 


Table 6 Deposition environments for barchan dunes at six sample sites in the Yarlung Zangbo River basin 


R Deposition environment (% of total particles) 
Sample site 


Aeolian Beach Shallow marine Turbidity current Fluvial 
S1 - - 10.3 7.7 2.6 
S2 7.7 - 5.1 - 7.7 
$3 - - - - 20.5 
S4 - - 10.3 5.1 2.6 
S5 - - - - 20.5 
Total 7.7 0 25.7 12.8 53.9 


Note: No data is available for S6. -, no data available. The grain size parameters of the dunes (Table 5) were used to determine the 
deposition environments based on the Sahu (1964) discriminant function. 
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As described above, the characteristics of the sediments were generally controlled by the wind 
environment, with a certain amount of spatial variation. However, due to the influence of 
topography and material sources, the effect of wind on the aeolian sands was more difficult to 
describe than in lower-elevation sand seas in northern China. 


6 Conclusions 


In this study, we modified the threshold wind velocity to account for the lower air pressure at our 
high-elevation study sites in the Yarlung Zangbo River basin of the southern Tibetan Plateau. We 
analyzed the spatial and temporal variation in the near-surface wind regimes using data from 11 
meteorological stations in and around the study area, and examined the relationships between this 
variation and the aeolian sediment environment. 

We found that the wind environment showed strong spatial and temporal differences. Wind 
speeds and the frequency of sand-driving wind tended to decrease from west to east, and were 
generally high in spring and winter, and generally highest in spring, but were generally lowest in 
summer and autumn. The mean wind speed has decreased significantly over time, but with regional 
and seasonal differences in the magnitude of the decrease. This variation of the mean wind speed 
was mostly caused by changes in the wind frequency. The maximum wind speed accounted for 
62.6% of the mean wind speed, which suggests that the maximum wind speed was responsible for 
weakening or strengthening the wind regime. DP and related parameters showed high spatial and 
temporal variation. DP and RDP generally decreased from west to east; RDP/DP was mostly 
intermediate (moderate variability in direction); and RDD mainly indicated winds towards the 
northeast. DP and RDP showed significant exponential decreases since 1980; in contrast, there 
were no significant temporal trends for RDP/DP and RDD. DP and RDP were generally higher in 
spring and winter, and were generally lower in summer and autumn; RDP/DP was generally lowest 
in winter, and was highest in summer, but the seasonal variation in RDD was complex. These 
variation characteristics indicated that the wind environment most likely responded to similar 
driving forces in the river's upper and middle reaches, but was different in the lower reaches. 

The wind regime corresponded well to the observed aeolian geomorphological patterns. DP 
represented the dynamic conditions for the development of the aeolian sediments, and RDD 
controlled the orientation of the aeolian sediments, but the sand drift direction was altered by 
variations in the local topography. The grain size of the sediments showed spatial variation, and 
was affected by the sand-driving wind speed. With increasing wind speed, the content of very fine 
and fine sands decreased, whereas the content of medium and coarse sands increased. However, 
the grain size characteristics of the aeolian sand showed strong influences of the local supply and 
relatively weak reworking by the wind. 
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